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Abstract 

Periodic boundary conditions are applied to a ferromagnetic spin lattice. A symmetrical lattice 
and its contributions all over space are being used. Results, for the Ising model with ferromagnetic 
interaction that decays as a l/r D+v law, are discussed in the mean field approximation. 
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Thermodynamics describes the behavior of systems with many degrees of freedom after they have 
reached a state of thermal equilibrium. However, their thermodynamic state can be specified in terms 
of a few parameters, they are called extensive or intensive state variables. At equilibrium, macroscopic 
observables are linear functions of iV (number of particles). If the function f(N) is an observable, the 
extensive property imposes f(N) = Nfw, where /at = f(N)/N and f(XN) = Xf(N) for very large 

Small size of system are considered for studying thermodynamic behavior with computer simula- 
tions, typically N ~ 500 particles. Available memory and speed of program execution (in the host 
computer) can limit the size of the simulated systems. However, in thermodynamics is necessary to 
obtain the trend of the thermodynamic properties for very large systems, iV — ► oo. Then, how can 
we expect to simulate with finite systems some properties related to infinite systems? The answer for 
this question depends on the property under investigation and the range of microscopic interparticle 
interactions. 

Periodic summations of forces do not depend on the nature of interactions, even for long range 
interactions, forces converge very fast [2]. If the range of interaction is smaller than the size of the 
system, computer simulations give the typical trends of the thermodynamic behavior of systems. When 
interactions are long-ranged, surface effects appear and become very important for all system sizes. A 
first approach requires to increase to very large size of the box. 

A long-range interaction is often defined as one in which the space interactions falls of no faster 
than l/r D , where D is the space dimension of the system. 

To eliminate surface effects from the computation let us use a particular procedure. Surface effects 
can be ignored for all system sizes if we use periodic boundary conditions. In periodic boundary 
conditions, a central cell is replicated throughout space to form an infinite lattice. In the course of 
the simulation, when a molecule moves in the central cell, its periodic image in every one of the other 
cells moves with exactly the same orientation in exactly the same way [2, 3, 4, 5]. Thus, as a molecule 
leaves the central cell, one of its images will enter through the opposite face. There are no walls at the 
boundary of the central cell, and the system has no surface. The central cell simply forms a convenient 
coordinate system for measuring locations of the N particles. 

Therefore, the size of the systems should depend on two parameters, M and M, as follows 

N = (2M + l) D Af (I) 
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where TV is the number of particles in every cell and (2M + 1) D represents the number of replications 
on all space. So, thermodynamic limit {N — > oo) is reached when M — > oo, but this solution is im- 
practicable even with modern computers, the time taken to evaluate the energy increases proportional 
to M 2 . An alternative way to approach the thermodynamic limit is to take a finite value of M and 
M -> oo. 

In recent years, due to the difficult to obtain exact solutions in systems with long-range interac- 
tions, much effort has been devoted to handling this kind of interactions in computational systems 
by molecular dynamics and Monte Carlo simulations. Another reason for developing simulation tech- 
niques is to test some approximated results [6]. In the present work, a long-range Ising lattice is 
studied in the mean field approximation. Furthermore, we look for a cutoff for involved summations 
in thermodynamics by a nonextensive scaling. 

We consider a computational cell in D dimensions (normally D=l,2,3) with size N and periodic 
boundary conditions given by repetition of a central cell to infinity. The central cell contains M spin^ 
Ising ferromagnets with two-particle potential[3, 7, 8] and magnetic moment fi that means, a system 
described by the Hamiltonian 

H = ~ E J (fij)SiSj -fihJ^Si (2) 
where h is an external magnetic field, Si = ±1 Vz and is the distance between two sites and 

M -. 

j(f) = j y — =, — (3) 

^ If + K\ D + U 

K=-M 1 1 

where J > measures the strength of the coupling. In precedent works, topics of systems with 
long-range interactions have been discussed, by several authors[6, 10]. A common behavior has been 
conjectured [9] for generic systems with arbitrary (long or short) range interactions. In this way, it is 
defined for continuos variations N* = {N~ V I D - l)/(-u/D) [9] and for a spin lattice N = 2 D+U N* [7]. 
Now, we can expect that quantities like internal energy, free energy, etc., per particle scale with N* 
and N, respectively . In the present case, variables include scaling with N; more explicitly, the so 
called intensive variables (T, P, etc.) scale with N Now, we calculate average of hamiltonian summing 
over the configuration space, the dependence on N 2 appears for an approximated evaluation of energy. 
Close to mean field approximation, the hamiltonian of a D-dimensional spin lattice with N lattice 
sites can be written 

N 

H = -J2E(N,h)S l (4) 
i=i 

where E(N, h) = —^N(S) — /j,h, the factor \ ensures that we do not count the same pair of spin twice. 
The quantity (S) is the average spin per site. This assumption helps us to recover all the results 
for system with long range microscopic interactions in a similar manner to the traditional mean field 
approximation. 

So, the average spin per site of the lattice is given by 



(S) = tanh 



i[^NJ{S)+tih 



(5) 



If h = 0, the magnetization Nfj,(S) will be zero for high temperature paramagnetic phase and it will be 
nonzero at lower temperatures where the spins have spontaneously aligned. This transition (critical) 
point is obtained when AJ/2fcsT = 1. So, the critical temperature closes to the mean field value 

T c = NJ/2k B (6) 

where T c /N takes a constant value, that does not depend on the range of the microscopic interaction 
in the present approximation. 
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In a simple numerical test for the critical point of a long-range Ising ferromagnet system , in one 
dimension, is obtained a departure from mean field approximation. We note that if v — > —1, the 
mean field model closes to the present system; but if v > — 1, simulated critical temperature T c /N is 
less than the mean field value. We have performed numerical calculations to test different ranges of 
parameters for system with long-range interactions and a more detailed study of critical phenomena 
for systems with this kind of microscopic interactions is in advance [11]. 

The internal energy is given by 



U(N) = ~NNJ{S) 2 . (7) 



Now, we can replace Eq.(l) in Eq.(7) 



U({2M + l) D Af) = (2M + 1) D+ H[/(A0, (8) 

where v < 0. So, U (N) violates the extensive property and it is a homogeneous function of degree 
(1 + H/D). 

The heat capacity becomes 

2Nk B (S} 2 (T c /T) 2 
cosh 2 ({S)T C /T)-T C /T 
and the magnetic susceptibility in the limit h — > is 



C = , — " ' ' (9) 



x = K 2 Jm (10) 

N J cosh 2 ( (S) T c /T) -T c /T V ; 

Finally, it is important to understand the microscopic behavior of the particles, because computation 
of the mesoscopic and macroscopic thermodynamical properties comes from the averages of quantities 
per particle when the number of particles increases to infinity. Advances in calculation of these typical 
and physical quantities is very important in the study of phase transitions and critical phenomena for 
this kind of systems. 
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